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In this paper an equation of state of neutron star matter which includes strange baryons in the 
framework of Zimanyi and Moszkowski (ZM) model has been obtained. We concentrate on the 
effects of the isospin dependence of the equation of state constructing for the appropriate choices 
of parameters the hyperons star model. Numerous neutron star models show that the appearance 
of hyperons is connected with the increasing density in neutron star interiors. Various studies have 
indicated that the inclusion of S meson mainly affects the symmetry energy and through this the 
chemical composition of a neutron star. As the effective nucleon mass contributes to hadron chemical 
- - - potentials it alters the chemical composition of the star. In the result the obtained model of the star 

^-S I not only excludes large population of hadrons but also does not reduce significantly lepton contents 

. in the star interior. 

■ PACS numbers: 24.10.Jv, 26.60.+C 
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\0 : I. INTRODUCTION 

Compact stars which are known from observations can be classified into two distinctive groups. The first one is 
^ ' exemplified by white dwarfs and the second by neutron stars. Neutron stars themselves are identify with pulsars 
' and compact X-ray sources At the core of a neutron star the density of matter ranges from a few times of the 
density of normal nuclear matter to about an order of magnitude higher. Thus various exotic forms of matter such as 
I hyperons or quark- hadron mix phase are expected to emerge in the interior of a neutron star The appearance of 
. these additional degrees of freedom and their impact not only on neutron stars but on proto-neutron stars structure 
' and evolution as well has been the subject of extensive studies |2|. Properties of matter at such extreme densities 
, are of particular importance in determining forms of equations of state relevant for neutron stars and in examining 
their global parameters |^ . Theoretical description of hadronic systems should be performed with the use of quantum 
chromodynamics (QCD) as it is the fundamental theory of strong interactions. However, at the hadronic energy scale 
where the observed experimentally degrees of freedom are not quarks but hadrons the direct description of nuclei 
O , in terms of QCD become inadequate. Other alternative approach has to be formulated one of which is quantum 
^ i hadrodynamics (QHD) giving quantitative description of the nuclear many body problem. QHD is a relativistic 
^ ' quantum field theory in which nuclear matter description in terms of baryons and mesons is provided. The original 
model (QHD-I) contains nucleons interacting through the exchange of simulating medium range attraction a meson 
and uj meson responsible for short range repulsion. Extension (QHD-II) of this theory |8[ |9[ |oj includes also the 
^ . , isovector meson p. Nonlinear terms in the scalar and vector fields were added in order to get the correct value of 
' the compressibility of nuclear matter and the proper density dependence in the vector self-energy. The variation of 
nucleon properties in nuclear medium is the key problem in nuclear physics. In order to incorporate quark degrees 
of freedom in the analysis of nuclear many-body system Guichon provides a quark- meson coupling model (QMC). 
The extension of the QMC theory namely the quark mean field model (QMF), describing a nucleon with the use 
of the constituent quark model, has been successfully applied to study the properties of both nuclear matter and 
finite nuclei. The model considered in this paper is an alternative version of the Walecka approach with enlarged 
meson sector. In the interior of neutron stars the density of matter could exceed normal nuclear matter density up 
to a few times, in such high density regime nucleon Fermi energies exceed the value of hyperon masses and thus 
the new hadronic degrees of freedom are expected to emerge. The higher the density the more various hadronic 
species are expected to populate. The onset of hyperon formation depends on the hyperon-nucleon and hyperon- 
hyperon interactions. Hyperons can be formed both in leptonic and baryonic processes. Several relevant strong 
interaction processes proceed and establish the hyperon population in the neutron star matter. Neutron star models 
are constructed at different levels of complexity starting from the most elementary one which assumes that neutrons 
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are the only component. The more sophisticated version is formulated under the assumption that the neutron star 
matter has to obey the constrains of charge neutrality and /? equilibrium. Thus the model considered describes high 
isospin asymmetric matter and it has to be extended by the inclusion of isovector-scalar meson ao (980) (J meson) [0. 
For the sake of completness additional nonlinear vector meson interactions are included. When strange hadrons are 
taken into account uncertainties which are present in the description of nuclear matter are intensified due to the 
incompleteness of the available experimental data. The standard approach does not reproduce the strongly attractive 
hyperon-hyperon interaction seen in double A hypernuclei. In order to construct a proper model which do include 
hyperons the effects of hyperon-hyperon interactions have to be taken into account. These interactions are simulated 
via (hidden) strange meson exchange: scalar meson /o(975) {a* meson) and vector meson 0(1020) {(f> meson) and 
influence the form of the equation of state and neutron stars properties. 

The solution of the presented model is gained with the mean field approximation in which meson fields are replaced 
by their expectation values. The parameters used are adjusted in the limiting density range around saturation density 
Pa and in this density range give very good description in finite nuclei. However, incorporation of this theory to higher 
density require an extrapolation which in turn leads to some uncertainties and suffers of several shortcomings. The 
standard TMl parameter set for high density range reveals an instability of neutron star matter which is connected 
with the appearance of negative nucleon effective mass due to the presence of hyperons. The Zimanyi-Moszkowski 
(ZM) fl^ model in which the Yukawa type interaction gsNf is replaced by the derivative one {gsN^ /MN)iJ'Njiyd'^4' 
exemplifies an alternative version of the Walecka model which improves the behaviour of the nucleon effective masses. 
It also influences the value of the incompressibility K of neutron star matter. The derivative coupling effectively 
introduces the density dependence of the scalar and vector coupling constants. Knowing the form of the equation 
of state (EOS) is the decisive factor in determining properties of neutron stars such as: central density, mass-radius 
relation, crust extent or the moment of inertia. 

The essential goal of this paper is to obtain within the described above model the equation of state for the neutron star 
matter on the basis of calculations carried out for asymmetric nuclear matter in the relativistic mean field approach 
(RMF) and to compare the obtained results with ones that are relevant for high density calculations - namely with 
a quark star. The inclusion of S meson affects the neutron stars chemical composition changing the proton fraction 
which in turn affects the cooling mechanism. If the proton fraction is higher than the critical value of about Yp ^ 0.11 
the direct URCA processes can proceed and this enhances the rate of neutron star cooling. Whether the proton 
fraction can exceed the critical value and at what density it occurs depends on the model. The proton fraction is 
almost entirely determined by the isospin-dependent part of the EOS thus the inclusion of d meson and nonlinear 
vector meson interactions infiuence the neutron star structure and properties. 

This paper is organized as follows. Section 2 outlines the extended model with derivative coupling including hyperons 
and additional mesons, together with the collected equations of motions. Their solutions enable the construction of 
the equation of state . In section 3 the equilibrium conditions leading to the relevant hyperon star composition are 
presented together with the chosen values of parameters. Section 4 contains numerical results and the discussion of 
their influence on neutron stars properties. 



For the description of properties of the infinite nuclear matter with nonzero strangeness the relevant parts of 
the SU(3) structure have been involved. Due to the parity conservation the appearance of pseudoscalar mesons is 
forbidden. The Lagrangian function for the system can be written as a sum of a baryonic part including the full octet 
of baryons together with baryon-meson interaction terms, a mesonic part containing additional interactions between 
mesons which mathematically express themselves as supplementary, nonlinear terms in the Lagrangian function, and 
a free leptonic part 



II. THE MODEL 



(1) 



The interacting baryons are described by the Lagrangian function Cbm which is given by 
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where the spinor = (V'Af, V'A, V'Si is composed of the following isomultiplets 



^'S = I '/'so I , *H = ( ° 



The covariant derivative I?^ is defined as 



Df^ = di_^ + igujBi^i^i + ig^Bi'i^i + igpBhBT°'b'^^- (3) 

The model considered represents an alternative version of the Walccka model in which the Yukawa interaction term is 
replaced by the derivative coupling one. This Zimanyi-Moszkowski ||l2| method allows to solve the problem of too low 
effective nucleon masses achieved in the original approach. Rescaling the baryon field in a way proposed by Zimanyi 
and Moszkowski JT^ the modified Lagrange function for interacting baryons is obtained 

r ^1 ■ t^r, I /'i , 9sB(y + ga'B(T* + hBgsBT°-5°-\'^ - 

i^B = -2^'>lJBiTD^ip ~ — 1 iPbMbiJjb- (4) 

B B ^ ' 

The mesonic part of the Lagrangian function (Q) is given by 

+ (gpg,)2A4&>"^^2 - -^R^R^^" + \mXb'^^ + \c,{u:,uo^f + \ahlb'^n^ . (5) 
The field tensors i?^^,, fi^^, are defined as 

R^^d^bl^d^bl+gpSabcblbl, (6) 



^^j,u ^ df^uju - duujf^, (j)f_,u ^ df_,(t)i, ~ d^(j)fj,. (7) 

The potential function U{a) possesses the well-known polynomial form introduced by Boguta and Bodmer jo) 

U {<J) = ^mla^ + igsa^ + jgia*. (8) 

The baryon mass is denoted by Mb whereas (i — a,u!, p, a* , (f), S) are masses assigned to the meson fields and they 
are taken at their physical values. The parameters entering the Lagrangian function are the coupling constants g^B, 
gpB, gsB, gsB, g<T*B, g4>B for meson fields and the self- interacting coupling constants ga, 34, C3, ^. These parameters 
are adjusted to reproduce the bulk properties of the symmetric nuclear matter at equilibrium. They are collected in 
Table 2. As the 5 meson field and nonlinear vector meson interactions carry isospin they contribute to the symmetry 
energy Es of the system. This imposes constrains on gp and 175 A4 A„ and ^ parameters. They are adjusted to get the 
experimental value of the symmetry energy Es which is equal 32 ± 4 MeV . Thus the parameter gp has been redefined 
by comparison with the one in the original TMf parameter set. 

The presence of hyperons demands additional coupling constants (they are collected in Table 3) which have been 
established with the use of experimental data and theoretical analysis. 

The inclusion of hyperons improve our understanding of neutron stars at higher densities. There is common belief that 
at suitable densities a quark star is more stable than the neutron one. Thus it is interesting to compare the properties 
of hyperon stars with those obtained in the framework of QMF model for a quark star |6|. Now the Lagrangian 
function possesses the following form 

Cqmf = q[i^^Dp, - nicjq + ^dpipad^ipa - ■^mlipafa (9) 
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where q denotes a quark field with three flavors, u, d and s, and three colors. The construction of the model mimics 
the relativistic mean field theory, where the scalar a and the vector meson lu fields do not couple with nucleons 
but directly with quarks. The quark mass has to change from its bare current mass due to the coupling to the 
a meson. In the framework of QMF model with the use of the TMl parameter set quarks constituent masses are 
rric^u = "rric^d = 367.61 MeV and mc,s — 504.1 McV whereas the bag constant is chosen at the level of S^/'* = 154.5 
MeV. There are no significant differences between the presented above Lagrangian function (|o|) and the Lagrangian 
function (|^). They differ in the baryonic sectors whereas the mesonic sectors for both models are nearly the same. 
There are only differences in coupling constants: ga-q = gaB/i,gi^q = gujB/'^,gpq = gpB,gsq = gsB- We restrict 
ourselves to the isospin SU(2) unbroken symmetric case, m" = mj?, so 




cjSfj' = rn^.d (11) 



Generators of the U(3) algebra A° = {A° = 1/2/8/, A*} (where / is an identity matrix, A'are Gell-Mann matrices 
SU(3) algebra) obeying Tr{y\^) = 25ab- Restricting only to U{2) x U{1) subalgebra (a = {0,1,2,3,8}) case the 
simplest version of the QMF theory can be obtained. Defining the new base with t", a = {0, 1, 2, 3, 4} as 



r" = I 1 I = ( 0' ) M * = {1>2,3}, r4 = x/2 I 1 (12) 



1 
1 


the meson fields may be decomposed as follows 



A /O 0\ 

for i = {1,2,3}, = \/2 

/ \ 1 y 



and 



Vp = VI^X^ = LO^T^ + , 



_4 



Here the QMF model is enlarged by the inclusion of the isovector 5 meson. Similarly to nucleons it splits 
u and d masses. Both 5 and p mesons may be neglected in the case of the symmetric nuclear matter. 

Baryon and meson parameters 



Particle 


m [McV] 


J I S 


Y 


h 


Q 


n 


939.566 


1/2 1/2 


1 


-1/2 





P 


938.272 


1/2 1/2 


1 


1/2 


1 


A 


1115.63 


1/2 -1 











E- 


1197.43 


1/2 1 -1 





-1 


-1 


YP 


1192.55 


1/2 1 -1 











E+ 


1189.37 


1/2 1 -1 





1 


1 


S° 


1314.9 


1/2 1/2 -2 


-1 


1/2 







1321.3 


1/2 1/2 -2 


-1 


-1/2 


-1 


cr 


550 
















783 


1 
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770 


1 1 





-1,0,1 


-1,0,1 


5 


980 


1 





-1,0,1 


-1,0,1 


a* 


975 
















1020 


1 












Table 2 

Parameter sets used in this paper 



Parameter 


TMl 


ZM 


ZM + 5 


ZM + () + nonlinear terms 


M 


938 McV 


938 McV 


938 McV 


938 MeV 


C3 


71.3075 











53 


7.23 /m-i 











9i 


0.6183 













12.6239 


6.671 


6.671 


6.671 


9aN 


10.0289 


7.8449 


7.8449 


7.8449 


9pN 


9.2644 


8.9 


9.5 


9.5 


9SN 








3.1 


3.1 


9<7'N 














94>N 














Av 











0.008 


A4 











0.001 


C 











0.5 



Table 3 
Hyperon-meson couplings 



Parameter 


Value 


Parameter 


Value 


9h<y 


0.5207 g^N 




0.5815 g^N 


9Ap 









V2 „ 


9Au> 


3 9u}N 


9A4> 


9T.a 


0.1565 QaN 


.9s<T* 





9Sp 


2 9pN 




2 gsN 

-— 9ojN 


9T.UJ 


3 9u}N 




9Ea 


0.2786 g^N 


9Ea' 


0.5815 g^N 


9Sp 
9-s.uj 


9pN 
^ 9ujN 


gss 

9S4, 


gsN 

2\/2 „ 
— 9uN 
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As it is usually assumed in quantum liadrodynamics the mean field approximation is adopted and for the ground 
state of homogeneous infinite matter quantum fields operators are replaced by their classical expectation values. Thus 
one can separated mesonic fields into classical mean field values and quantum fluctuations which are not included in 
the ground state: 

(13) 





= a - 


\- S 


a* 


= a* 


+ s* 


6" 


= r 


+ (15^" 






+ foS^o 
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+ WO^nO 




= K 





■3a 



The derivative terms are neglected and if one assume homogenous and isotropic inflnite matter only time-like com- 
ponents of the vector mesons will survive. The fleld equations derived from the Lagrange function at the mean fleld 
level are 

mis + gzs^ + qas^ - 2{gpBgsBf ^irls = ^gsBM^j:f ,,^S{Meff,B,kF,B) (14) 

B 

mlwo + csivl + 2{gpBgujBfKrlwo = ^ g^BriB (15) 

B 

mlro + 2{gpBgsBf^Aros'^ + "^{gpBQuBfKrowl = ^ gpBhBriB (16) 

B 

mjd^ = J2 96Bl3BSiMeff,B,kF,B) (17) 
B 

=Y,9a'BMlff^^^S{M,ff^B, kF,B) (18) 

B 

mlfo^J^d'PBnB. (19) 

B 

The function S{Meff^B, kF,B) is expressed with the use of an integral 



(20) 



where Jb and Izb are the spin and isospin projection of baryon B, kF,B is the Fermi momentum of species B, hb 
denote the baryon number density. The Dirac equation for baryons that is obtained from the Lagrangian function 
has the following form 

[i^^dp, - MB,eff - 5a;B7% - 5<^b7°/o - ^5pB7°^^r-o) V = (21) 

with MB,eff being the effective nucleon mass generated by the nucleon and scalar fields interactions and is defined 
as 

_ ^B ^22) 

B,eff - ^ _^ ^^^^^ _^ ^ hBg5Bd)/MB ^ ' 
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III. THE EQUILIBRIUM CONDITIONS AND COMPOSITION OF MATTER 

In the high density regime in neutron star interiors when the Fermi energy of nucleons exceeds the hyperon masses 
additional hadronic states are produced. The onset of hyperon formation depends on the attractive hyperon-nucleon 
interaction. The higher the density the more various hadronic species are expected to populate. They can be formed 
both in leptonic and baryonic processes. In the last one the strong interaction process such as 

n + n ^ n + A (23) 

proceeds. There are other relevant strong reactions that establish the hadron population in neutron star matter e.g.: 

A + n^j:^+p A + A^S^+p (24) 

The comparison of weak interaction time scales (10"^" s) and the time scale connected with the lifetime of a relevant 
star indicate that there is a difference between the matter in high energy collisions which is constrained by the isospin 
symmetry and the strangeness conservation whereas neutron star matter by the charge neutrality and the generalized 
/3-equilibrium. Thus realistic neutron star models describe electrically neutral high density matter with no strangeness 
conservation being in /3 equilibrium. The last condition implies the presence of leptons. Mathematically it is expressed 
by adding the Lagrangian of free leptons to the lagrangian function (Q) 

5] V^i(«7^a^-m,)V/- (25) 

Neutrinos are neglected here since they leak out from a neutron star, whose energy diminishes at the same time. After 
electron chemical potential fi^ reaches the value equal to the muon mass, muons start to appear. Equilibrium with 
respect to the reaction 

e~ i-^ fi^ + Ve + (26) 

is assured when /i^ = /Xe (setting fi^^ = fj,p^ = 0). The appearance of muons reduces number of electrons and affects 
also the proton fraction. The introduction of the asymmetry parameter fa which describes the relative neutron excess 
defined as 

fa = ^^^^ (27) 

71 AT 

allows to study the symmetry properties of the system. The equilibrium conditions between baryonic and leptonic 
species which are present in the neutron star matter lead to the following relations between their chemical potentials 
and constrain the species fraction in the star interior 

^J■p = A*s+ = ^J■n - ^J■e = a*s" = = Mn (28) 

where the baryonic chemical potential is given by the relation 



fJ'i^ y kp B + ^B,ef f + 9^3(^0 + g<t>Bfo + hsgpETo (29) 

Similarly to the asymmetry parameter fa a parameter which specify the strangeness content in the system and is 
strictly connected with the appearance of particular hyperon species in the model has been introduced. 

= + + . (30) 

The RMF theory as an effective one requires the knowledge of coupling constants. There are several selected param- 
eterizations of the theory strictly connected with the assumption that the exchange of scalar, isoscalar-vector, vector 
mesons and two hidden-strangeness scalar and vector mesons are responsible for interactions between specific con- 
stituents of the matter. Particular sets of coupling constants are associated with the description of nucleon-nucleon, 
hyperon-nucleon and hyperon-hyperon interactions and are indispensable for the the construction of the equation of 
state which in turn is applied for determining a neutron star properties. The enlarged Walecka model uses the TMl 
parameter set (Table 2)p^. More realistic description of a neutron star requires taking into consideration not only 
the interior region of a neutron star but also remaining layers, namely the inner and outer crust and surface layers. 
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The composite EOS constructed by adding Bonn [Q and Negele-Vautherin (NV) equations of state (describing 
the inner crust) to the TMl one allows to calculate the neutron star structure for the entire neutron star density 
span. The parameters describing the nucleon-nucleon interactions are created in order to reproduced the properties of 
the symmetric nuclear matter at saturation such as the binding energy, symmetry energy and the incompressibility. 
Other parameters are related to the hyperon-nucleon interactions. The scalar meson coupling to hyperons can be 
calculated from the potential depth of the hyperon in the saturated nuclear matter 

up (po) = -gsVP + 9u,Y^- (31) 



The vector coupling constants for hyperons are determined from 5C/(3) symmetry as |19 



11 1 

-j^g^A = 2^^^^ ^ ^ gSt^Af (32) 
1 

1^9 py. = 9pB. = gpN] gpA = 

o _ _ 2%/2 
2g</>A — — ,g0H — —^g,^N- 

The single A potential in nuclear matter is well determined as (||2^) —Uj^{po) = 27 — 28 MeV. Recent analysis of 
atomic data ||2^, ^ indicate for repulsive S potential in the interior of nuclei 1/^ (po) = 30 MeV thus S hyperons 
do not appear in the bulk matter calculations. The interpretation of the S hyperons data gives the value of the 
potential well depth —U^ (po) = 18 MeV. The experimental data concerning hyperon-hyperon |jl^, ^ interactions 
are extremely scare. Analysis of events which can be interpreted as the creation of AA hypernuclei allows to determine 
the well depths of hyperon in hyperon matter 

-Uhipo) = 40MeF (33) 



Having specified parameters of the model the equation of state can be calculated. This has been done with the use of 
the energy-momentum tensor T^^ defined as 

T,.^2^-g,X. (34) 

It allows to calculate the pressure P and energy density e of the system. The pressure P is related to the statistical 
average of the trace of the spatial component of the energy-momentum tensor P = ^ < Tu >, whereas the energy 
density e equals < Too >■ Thus the complete form of the equation of state includes contributions coming from meson, 
fermion and baryon fields and finally one can get the following equations for the energy density e and pressure P 



2 " 4 
with eb and sl given by 



-TOp(ro)^ + ^rns(f + -m^f^ + -m^. {s*f (35) 

1 2 2 I 3 4 , 222 2|A 2222|1a44i , 

™c^^"o + + U[s) + SA-^g^g^roWg -I- kig^g^r^Lp + -CSp^-o + (b + 



^5 = E i k'dk^/^TM!^ (36) 



p 



(37) 
(38) 



im,. {s*f - u{s) + Kglglrl^l + t.^glglrlv^ + ^C5>o + Pb + Pl 
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k'^dk 



(40) 



IV. NEUTRON STAR PARAMETERS 



Besides other properties of neutron stars the value of their masses and radii are very sensitive to the chosen model of 
strong interactions which in turn lead to the significant constrains on the form of the equation of state of the neutron 
star matter. In this paper three different groups of parameters have been applied in order to examined neutron star 
properties. In all of them the nucleon-hyperon E interaction is assumed to be repulsive. The first case marked as set I 
does not contain 5 meson, the second one (set II) do include S meson, in the third both 5 meson and nonlinear vector 
meson interactions are taken into account (set III). As the neutron star matter is of sizable asymmetry the last case 
seems to be the most adequate for the complete description of the asymmetric neutron star matter. The main effect 
of such an extension of the theory becomes evident studying properties of the neutron star matter especially baryon 
masses splitting and the form of the equation of state. Fig.|l| depicts the effective baryon masses obtained for the 
three mentioned above cases as functions of the baryon number density ub- For the second and third cases there are 
noticeable in-medium baryon masses splitting for each isomultiplet. This effect reaches the maximum value for the 
third group of parameters, for the density span of (3 — 4)yOo and depends on the considered baryonic masses. Depending 
on the sign of the third component of particular baryon isospin the S meson interaction increases the proton and 
effective masses and decreases masses of neutron, and S^. Contrary to this situation, for parameters marked as 
the set I, the baryonic masses for the given isomultiplet remain degenerate as it is shown in the left panel of Fig.|l](a). 
Having obtained the effective baryonic masses one can compare them with those obtained in the Walecka model. In 
the original Walecka approach the nucleon effective mass rapidly diminishes its value passing through zero and even 
become negative for higher densities. Throughout the effective baryonic masses the asymmetry of the system alters 
baryon chemical potentials what realizes in characteristic modification of the appearance, abundance and distributions 
of the individual fiavors. This is evident comparing the results obtained for the three mentioned above parameter 
sets. Fig. ^ presents fractions of particular baryon species Yb as functions of baryon number density Ub- Starting 
the analysis of these graphs from low and moderate densities it is evident that at very low density neutrons are the 
most abundant baryon and a star resembles the pure neutron one. For higher densities protons and electrons and 
then muons emerge. Fig. 2 points that the first strange baryon that appears at p = 2.5 x po is the A, it is followed by 
S" and S". In the presence of 6 meson and in the case when nonlinear vector meson interactions are included the 
sequence of appearance of hyperons is the same as in the first case, however shifted towards higher densities. Due to 
the repulsive potential of E hyperons their onset points are possible at very high densities which are not relevant for 
neutron stars. Taking into account results obtained when the attractive nucleon-hyperon S potential is assumed 
appear as the first strange baryon and than A and S+. emerge at the density above 7po, whereas other hyperons 
at even higher densities. 

The appearance of charged hyperons permits the lower lepton contents and thus charge neutrality tends to be guar- 
anteed without lepton contribution. This kind of deleptonization in the case of the set II and III parameterizations 
similarly to the baryon distributions takes place at higher densities. Larger effective baryon masses cause the shift of 
the given hyperon onset point especially for the charged ones towards higher densities. The emergence of H" hyperons 
through the condition of charge neutrality affects the lepton fraction and causes a drop in their contents. The proton 
fraction Yp has a crucial role in neutron stars cooling history. At a certain critical value of Yp the direct URCA 
processes for neutrino emissions are allowed. The condition which has to be satisfied for the URCA process is given 
by the relation between Fermi momenta of nucleons and electrons 



The threshold proton fraction for this process is about 0.11 in the case when the neutron star matter consists only of 
nucleons and electrons. Inclusion of muons increases K, to the value ~ 0.14. The proton fractions which are obtained 
for the presented parameter groups are shown in FigJS. The third parameter set gives the lower value of Yp. 



kpp -f kpf. > fci?; 



'n 



(41) 
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FIG. 1: Effective baryon masses as a function of baryon number density ub- Fig. la shows effective baryon masses obtained 
in the model without 5 meson (the parameter set I). In Fig. lb and Fig.lc the effective baryon masses are splitted due to the 
presence of 5 meson. This effect is even more significant in the case of more asymmetric matter (set 111). 
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FIG. 2: The equilibrium compositions as functions of baryon number density ub for given parameter sots. Fig. 2a shows rcsuhs 
obtained for the parameter set I. Fig.2a and 2b depict equiUbrium compositions of matter which contain 5 meson (set II) and 
additional vector meson interactions (set III). 
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FIG. 3: Proton fractions Yp, as functions of baryon number density ub- Solid lino represents the proton fraction for the most 
asymmetric matter. Dotted and dashed lines are obtained for parameter set I and II respectively. 
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Baryon distributions are strictly connected with the behavior of meson fields. The values of scalar and vector meson 
fields influence the onset points of individual hyperon species. Fig. ^ shows meson fields as functions of the baryon 
number density. 

The calculated form of the equation of state (EOS) determines the physical state and composition of matter at high 
densities and is presented in Fig.^j. The relative hadron-lepton composition calculated for all parameter groups can be 
analyzed through the density dependence of the asymmetry parameter fa and the parameter fs which is connected 
with the strangeness contents. Fig.|^ presents both parameters as functions of the baryon number density ub- As 
the density increases the asymmetry parameter decreases. The third parameter group gives the higher value of the 
asymmetry parameter. For parameter sets I and II the asymmetry is comparable. The strangeness contents increases 
with the density. The inclusion of 6 meson and nonlinear vector meson interactions results in the lowest hyperon 
contents. 

The obtained form of the equation of state serves as an input to the Oppenheimer-Volkoff equations and determines 
the structure of spherically symmetric stars. 

dPjr) Gmir)pir) (1 + gg)(l + ^) 

dr r2 1 _ 2Gm(r) ' ^*^) 

r 

dm(r) 2 / \ 

dr 



Numerical solutions of these equations allow to construct the mass-radius relations. These relations for the chosen 
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FIG. 5: The equation of state (EOS). Straight, dashed and dotted lines represent models constructed with the use of parameter 
sets 1,11 and 111 respectively. 
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FIG. 6: The asymmetry parameter /o and the strangeness parameter /s as functions of baryon number density ns. 
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FIG. 7: Mass radius relations. Asterisks represent points where direct URCA processes start whereas arrows the points where 
hiperons emerge. 



form of the equations of state are presented in Fig|^. Arrows represent points at which hyperons emerge whereas 
asterisks the configurations for which URCA processes start to proceed. For the parameter set III these two marks are 
in the same position. This figure also shows that 6 meson itself and the additional nonlinear vector meson interactions 
change the maximal mass but not in a significant way. However, there are visible differences in the maximal radius 
configurations. For comparison, the mass-radius relation for a quark star is enclosed. In Figjl the obtained neutron 
star masses as functions of central densities are presented. 
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FIG. 8: The neutron star masses as functions of the central density of neutron stars. 
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FIG. 9; Gravitational binding energy as a function of central density of the neutron star. 



White dwarfs and neutron stars are purely gravitationally bound compact stars. The gravitational binding energy 
of a relativistic star is defined as a difference between its gravitational and baryon masses. 

Eb,g = (Mp - m{R))^ (43) 

where 

M„ = 4^ r drr^l ~ ^S^)-^p(r) (44) 
Jo c^r 

Of considerable relevance is the numerical solution of the above equation for the selected EOS. Results are shown 
in Fig. |[ For the sake of completeness calculations have been done not only for the presented parameter sets 
but also for the original TMl parameter set and for a quark star. From this figure it is evident that at moderate 
densities the standard neutron star model calculated with the use of the TMl parameter set is the most energetically 
favorable configuration. At higher densities this standard RMF approach become unsatisfactory. There is room for 
the improvement of this theory by the introduction of other degrees of freedom. The comparison of the models which 
contain hyperons with the Walecka one is performed. At high densities the models with hyperons are more favorable 
than the model calculated with the use of the TMl parameter set. This is directly connected with the value of effective 
baryon masses, obtained with the use of the TMl parameters, which diminish considerably their values which after 
passing through zero can take even the negative value. This can be interpreted as a break down of the theory and 
has its confirmation in the form of the gravitational binding energy. Fig.|9| depicts the gravitational binding energies 
as functions of the density. From this figure it is evident that the models without hyperons are energetically favorable 
at the density range relevant for neutron stars. The situation changes for higher densities and now models with 
hyperons become energetically favorable. There arc theoretical suggestions about the existence of a quark star but 
Fig.^ exhibits that a quark star configuration should appear at very high densities. 

The gained solutions of the structure equations allows as to carry out a similar analysis of the onset point, abundance 
and distributions of the individual hadron and lepton species but now as functions of the star radius R. Comparing 
results obtained for the three presented above cases (parameter set I, II and III) one can come to the conclusion that 
in all the cases the assumption of the repulsive E interaction shifts the onset point of E hyperons to very high densities 
and they do not appear in neutron star interiors calculated in this model. Two characteristic configurations have 
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FIG. 10: The equilibrium compositions for the maximum mass configuration as a function of the star radius. Fig. 10a is 
constructed for the parameter set I, Fig. 10b for the parameter set II. Panel (c) represents results obtained for the parameter 
set III. 

been considered. Namely the one connected with the maximum mass configuration and the second with that of the 
maximum radius. The very compact hyperon core which emerges in the interior of the maximum mass configuration 
consist of 'EP , S~ and A hyperons only for the second parameter set. The hyperon population is reduced to A and 
S~ for the parameter sets I and III. This can be seen in Fig.^ and confirm by the behavior of meson fields in the 
neutron star. Especially by the hidden strange a* and (j) mesons which appear in the vicinity of the neutron star center 
(Fig.|l^ and Fig.|l^) . For the maximum radius configuration for aU groups of parameters hyperons do not emerge in 
the interior of neutron stars. 
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FIG. 11: The equilibrium compositions for the maxirrmm radius configuration as a function of the star radius. Panel (a) is for 
the parameter set 1 whereas (b) and (c) for parameter sets 11 and 111. 
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FIG. 12: The behavior of meson fields in the configuration with maximum mass as functions of the star radius. Panels (a), (b) 
and (c) are for parameter sets I, II and III respectively. 
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FIG. 13: The meson fields in the neutron star with maximum radius as a fijnction of star radius. The panel (a) is for parameter 
set 1, the panel (b) for parameter set II the panel (c) for parameter set III. 

V. SUMMARY AND CONCLUSIONS 

In this paper the complete form of the equation of state of hyperon matter has been obtained with the use of the 
derivative coupUng model in the framework of an extended RMF theory which besides hyperons and leptons includes 
the extended meson sector with additional 5 meson and hidden strange mesons a* and (j). The model considered 
is also supplemented with nonlinear vector meson interactions. This enlargement alters the symmetry properties of 
neutron star matter and through this neutron stars parameters. The value of baryon effective masses depend on the 
scalar meson condensates and at high densities when hyperon species appear the possibility of negative nucleon masses 
emerges. The derivative coupling model allows to avoid this difficulty reproducing reasonable value of baryon effective 
masses for densities relevant for neutron stars. The inclusion of 6 meson and nonlinear vector meson interactions 
influences the chemical composition of a neutron star. This is especially evident comparing the effective baryon 
masses in the density span (3 — 4) x po and equilibrium compositions of the star. For the third group of parameters 
the higher value of asymmetry has been obtained. This changes the properties of a neutron star diminishing the 
hyperon core extent. The asymmetry of the system also influence the star radius, for the third group of parameters 
the one can obtain the lower value of the star radius. There is also possible a configuration representing a star 
with lower densities (a maximum radius configuration) which excluded the existence of a hyperon core. The model 
considered excludes a large hyperon fraction which can be connected with thermal properties of a hot neutron star. 
As the most populated strange baryon is the A hyperon it does not reduce significantly number of leptons in the star 
interior and thus the models calculated with the use of the parameter set I and II do not exclude rapid cooling rate of 
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the star. This is even more evident analyzing the proton fractions obtained for all parameter groups. The equilibrium 
proton fraction is also determined by the nuclear symmetry energy. The parameter sets I and II permit higher values 
of proton fraction which is indispensable for URCA processes to proceed. However, the equilibrium proton fraction 
Yp is significantly reduced for the third group of parameter. 

Analyzing the gravitational binding energy one can come to the conclusion that configurations with hyperons are 
energetically favorable than the one obtained with the use of TMl parameter set for higher densities. 
All assumption which have been made namely: the derivative coupling model being connected with the higher effective 
baryon masses, the inclusion of 6 meson and nonlinear vector meson interactions, and the repulsive nucleon-hyperon 
S interaction lead to the neutron star model with the value of maximum mass close to 1.5 Mq with the reduced value 
of proton fraction and very compact hyperon core. The calculation of the quark matter equation of state, allows to 
construct the mass-radius relation for the quark star. Comparing gravitational binding energies one can come to the 
conclusion that the addition of hyperons to the model shifts the stable hyperon matter configuration towards higher 
densities even to the density range which is relevant for a quark star and at the same time makes the existence of a 
pure quark star more problematic. 
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